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Abstract: We consider the problem of sliding motion of a charge- density- wave subject 
to static disorder within an elastic medium model. Starting with a field-theoretical 
formulation, which allows exact disorder averaging, we propose a self-consistent ap- 
proximation scheme to obtain results beyond the standard large-velocity expansion. 
Explicit calculations are carried out in three spatial dimensions. For the conductivity, 
we find a strong-coupling regime at electrical fields just above the pinning threshold. 
Phase and velocity correlation functions scale differently from the high-field regime, 
and static phase correlations converge to the pinned-phase behaviour. The sliding 
charge-density-wave is accompanied by narrow-band noise. 
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1 Introduction 



The peculiar dynamics of charge-density- wave (CDW) systems has attracted attention 
for more than two decades |l|] . Among the most prominent transport properties observed 
is the non-hnear electric current response to a constant electric field, the pinning of the 
CDW below finite threshold, harmonic and subharmonic mode-locking in an additional 
oscillating drive, and hysteretic effects. It is well established that this behaviour has 
its very origin in the interaction of the CDW with quenched disorder in the parent 
substance. 

The CDW state is formed through the instability of the quasi one-dimensional electron 
system with respect to static density modulation below the Peierls transition, induced 
by electron-phonon interaction. It is characterized by a complex order parameter, the 
amplitude being proportional to the amplitude of uniaxial electron density modulation, 
and the phase modes as low-energy collective excitations describing deformations 01 . 
Microscopic theory [0] shows the classical character of phase dynamics, which is usually 
overdamped. Therefore, the intriguing difficulties of theoretical approaches to CDW 
dynamics stem from the many degrees of freedom coupled non-linearly to disorder. We 
limit ourselves to phase dynamics throughout this paper, and neglect a possible influence 
of topological defects as well. 

In the following we focus on the sliding motion of the CDW driven by a constant 
field. The simplest picture describes the dynamics by a single variable. Its equation of 
motion is assumed to be that of an overdamped particle in a one- dimensional periodic 
potential supplementary pulled by a constant force. (Note the analogy to the dynamics 
of Josephson junctions.) Though this model renders correctly qualitative features as 
the non-linear force-velocity dependence and even pinning, the periodic potential is an 
unknown input and a sinusodial potential (independent of particle velocity) yields quan- 
titavely incorrect results, especially for exponents. Also it does not give any information 
on correlation functions. 

A perturbative approach to the problem, including the spatial degrees of freedom of 
the phase, was performed by Sneddon et al. 0, inspired by earlier work on vortex lattices 
in type-II superconductors 0. The disorder, which is considered as perturbation to the 
overdamped motion, is modelled by a static random potential with vanishing mean. In 
three space dimensions, the authors found in second order a correction to Ohm's law 
proportional to the square root of the field. Generalizing to arbitrary dimensions (i, the 
perturbation expansion breaks down at small fields for d ^ A. 

The more recent view of the depinning transition as a dynamic critical phenomenon 
requires a dynamic field theory, i.e. a formulation, where disorder averaged dynamics is 
mapped exactly on a self- interacting system ^ 0. Thus the apparatus of renomalization 
group theory can be deployed §]. Interestingly, non-standard elements remain, e.g. the 
effective periodic force correlator has to be given assertively a cusp at its maximum. 

Recent numerical 0] and experimental work seems to reveal a more complex 
scenario above threshold. The aim of this paper is to formulate and work out self- 
consistent approximations within the dynamic formalism, which has several advantages 
compared to the standard approach []T^ 
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In the next Section we introduce elementary notions of phase dynamics and convenient 
reduced variables. Our starting point for further investigation of non-linear dynamics 
is the formulation of Eckern et al. [0, which is presented in Sec. 3. This functional 
integral representation of the generating functional combines the dynamic variables of 
forward and backward time evolution in independent external fields. Due to a pecu- 
liar normalization property, exact disorder averaging is feasible, resulting in an effective 
interaction. In Sec. 4 we prove a functional differential equation for the generating func- 
tional of connected Green's functions, on which self-consistent approximation schemes 
can be based. It is neither restricted to a particular interaction term nor to a definite 
field-theoretical formulation. We define especially the mean field and the self-consistent 
Hartree approximation. The latter is used in Sec. 5 to extract the behaviour of the CDW 
close to threshold. Finally conclusions are given in Sec. 6. In an appendix we mention 
a different line to systematic approximations based on functional Legendre transforms, 
applied e.g. in [jl4[, and prove equivalence of a specific approximation to one of ours. 



2 Phase dynamics 

The CDW represents a typical low temperature state of highly anisotropic materials. Due 
to the nesting property, the electron-phonon interaction induces correlated scattering of 
electrons between two opposite parts of the Fermi surface, leading to both a spatial 
modulation of the electron density and a lattice distortion oc cos{Qx) with Q = 2kp. 
The CDW direction is chosen along the x-axis, and kp denotes the Fermi wave number. 
Thereby a gap A for quasi-particle excitations is opened. In case of incommensurability 
with lattice vectors, the CDW is free to move along the x-axis. This can be described 
by a continous phase shift in the density modulation oc cos{Qx + ({)). Consequently, 
the CDW current density is related to the phase via j = enodt(p/Q, where uq is the 
condensate density. 

Omitting for the moment damping effects, the dynamics of the collective phase mode, 
allowing for slow temporal and spatial variations, is covered by the phase Hamiltonian 
or, alternatively, by the action 0, |l^ 

The first part includes inertia, elastic energies, and the electric field component E^. 
The condensate density calculated from microscopic theory reads no = HvfNfQ^I — Y), 
where vp, Np, and Y denote Fermi velocity, density of states per spin at the Fermi 
surface, and the Yoshida function. The prefactor of the time derivative is enhanced by 
the Frohlich mass ratio mp = 1 + 4A^/[A^^co'q(1 — Y)]. It depends on the dimensionless 
electron-phonon strength A, and the phonon frequency ujq. Typically, the transverse 
velocity v± is much smaller than vp. The second part, 

Su = nQ j dtd\Re[eexp(i^)], (2) 

contains the influence of disorder, which breaks translational invariance. The density 
uq = 2NpA/X represents the amplitude of the electron density modulation. In the 
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appropiate weak-pinning limit, disorder is simulated through static complex random 
fields C,- They can be chosen to have zero mean value, and correlations {C,C,) = { ^*^*) = 
and 

where tq is the characteristic time for scattering processes between opposite parts of 
the Fermi surface. Regarding damping effects, investigations of phase dynamics can be 
based, motivated phenomenologically, on the equation of motion 

no.a,.^^. (4) 

thereby introducing the damping constant 7. 

Typically, disorder leads to suppression of long-range order. This can easily be demon- 
strated in the present case by calculating the equal-time phase-phase correlations up to 
second order in ^. In three dimensions and specializing to the static case E^^ = 0, they 
are given by 

Mr, 0) - ^(0, 0)]^) = Vix/Lor + {rjL^y. (5) 

The quantities Lq = 4:'n-'^TQhNFVFV±nl/Q'^nQ and L± = \/v±Jvf Lq determine the length 
scale, where disorder becomes relevant. Apart from a numerical factor, this scale follows 
as well from minimizing elastic and disorder energy [0, |T^. The linear dependence on 
spatial separation persists to all orders in perturbation theory fl^. It has, however, 
become clear that perturbation theory must fail, because the pinned CDW can access 
a multitude of metastable configurations. More sophisticated approaches indicate a 
logarithmic growth with distance on long scales [|^, |20|, ^ . 

The quantities Lq, L±, to = 2QLI'j/vf, Eq = hvp/^eLl, and ^0 = JT'VFno/2QLlnQ, 
suitable for overdamped motion, are used to introduce reduced variables for the x- 
coordinate, transverse coordinates, time, electric field, and random fields. From now 
on we use exclusively these reduced variables, but maintain notation. For later con- 
venience some of the above equations are displayed again. The action is then written 
as 

f = a/ dtdV{^ id.^y - i (V^)2 + f/(v.) + E^^}, (6) 

where 

[/(y.) = Re[eexpM]. (7) 

Correlations are given by 

{ar)C{r'))=8n6{r-r'), (8) 
and the equation of motion reads 

fid^^ + dt^-V^^ = E, + U\^). (9) 
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The prefactor /i = mi?/2Q^i^o7^ typically small compared to one. Whereas A = 
LoL^no7 being large, indicates a strong suppression of thermal phase fluctuations. For 
numerical values, especially concerning NbSea, we refer to 0, |15|. 



3 Generating functional 



The essential building block of dynamical field theory [^, ^ is the contour along which 
time evolution of the system is considered: It runs from t = — oo to t = +00, and 
back to t = —00. For convenience, the contour integration is split into two ordinary 
time integrations by introducing an additional index a denoting forward {a = 1) and 
backward {a = 2) path. The phase variables are then combined into if"' = (y?^,^?^), 
where a = {a,t,r), as well as the external fields into i]"- = {if-,rf). Minus signs from 
backward integration can be hidden in lower index quantities, e.g. t]^ = {1]^, —rf). This 
is supplemented by the convention that only indices appearing as both upper and lower 
ones are subject to path summation, and time and space integration respectively. 

Taking an average with respect to the initial state, the complete dynamical informa- 
tion is included in the generating functional 

Zu[ri\ = J DM exp(i5o[¥^] + ^M^] + ^Va^')- (10) 

It is obvious that the exponent essentially contains the difference between the action (P) 
with variables of path one and two. Specifically, the action of the free system is 

So = ^^'^rV- (11) 

The free system is assumed to have been in thermal equilibrium in the infinite past, and 
the disorder contribution 

Su = Aj dtd^r{U{ifi^) - U{ip^)} (12) 

is switched on adiabatically. By construction, the generating functional is normalized to 
one for r]^ = yf, independent of Su. This most important feature can also be demon- 
strated perturbatively. 

A more transparent formulation is achieved by transforming to sum and difference 
variables ip"" via ip^ = ip°'Qj', the matrix Qj' being defined through 

= ^2 ^ ^i/2A, = - ^V2A. (13) 

The inclusion of A in the transformation takes into account the smallness of fiuctuations 
in if^ — (f"^. In order to keep covariance of the above expressions, we introduce new 
quantities fj^ = QaVb = QaQif^^cdi ^"^^ finally suppress from now on the 

overlines. Thus, using (0), the disorder term reads 



Su = 2K j dtdV sin((^V2A) U'i^"^) 



(14) 
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Omitting for the moment Su in (|T0|) , the functional integral is gaussian and can be 
calculated easily: 



ZoW = exp[(-i/2)r/,GS\]. (15) 

Both operators F^^, and Gq^ are related by r^fjCp'^ = 6^'^, with 5^"^ the identity. They have 
the matrix structure 



The normalization property of the generating functional, now fulfilled for 

Va=\ 1, (17) 



is a result related to the identity G^^ = 0. This generalizes to all n-point functions in 
the interacting case. The other entries, their Fourier transforms being displayed below, 
are the retarded and advanced components 



G^ = (/io;^ + i^ - qT\ = (G^y. (18) 

Finally the Keldysh component is given by the fluctuation-dissipation theorem (classical 
limit) 

G^ = {G^-Gt){e/u;), (19) 

where 6 = kBTto/Ah. In this formulation, the damping could be introduced easily. Also 
the separation of response and correlation functions and their possible independence in 
non-equilibrium is clearly visible. Thermal fluctuations of the phase have white noise 
character with amplitude F^ = 216'. In explicit calculations, in addition, we have to 
introduce a high-frequency cut-off in order to obtain finite local phase fluctuations. This 
is physically justified by the fact that the above action is valid only for low energy 
phenomena. 

The phase dynamics at fixed disorder potential U is characterized by the n-point func- 
tions {(f"''^ ■■ Given the generating functional Zjj, they follow as n-time derivatives 
at = 0, normalized to the functional itself at the same point. Due to the normalization 
property Zij[ri2 = 0] = 1, we may readily use the disorder averaged functional 

Z = {Zu) (20) 

to generate disorder averaged n-point functions y?"")) = {{ip^^ .. ip"-")^) . The av- 

eraging procedure can be performed inside the functional integral. Using the formula 
(exp(i5c/)) = exp(— (5(75(7 )/2) = exp(i5i), where the last equation defines the effective 
interaction 

Si = 2ni j dVdrdf {2Kf fAn[<p^ /2K) ?,m{>/ /2K) cos(^2 - i/), (21) 
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we finally obtain 



Z[ri]= / D[ip] exp(i5o[y?] + iSi[(p] + iri^ip''). 



[22] 



We remark that the field-theoretic approach to classical dynamics |]2^ , p5| , starting from 
(H), leads to the contributions ([I^) and ( ^1]) in the limit A ^ oo. 



4 Self-consistent approximations 

In the last chapter we derived an effective action that governs the disorder averaged 
phase dynamics. However it contains the non-linear term (PT|), and therefore one has to 
resort to an approximative analysis of the problem. The first step in our strategy is to 
write f^Tf) in the form 



5i = -iTr^exp(iCV). (23) 
Here the explicit structure is encoded in the symbols 

Tr^ = / d'r^dTAdfA {2tAeA)-\ (24) 



eA,eA = ±l/2A 



Kta - Ta) I ? - Kta - fA] ' 



1 ^ " "Ml 



(25) 



where A = (r^, ta, ta, e^, ca) combines internal variables being summed and integrated 
respectively. Then recalling that the integral of a derivative vanishes identically, we 
perform the derivative in /D[(y9] {S/Sip)^exp{iS[(p] + irjf^ip^) = to get 

J DM {{6/6ip),S[cp] + Va) exp(i5M + ivi>^') = 0. (26) 

A derivative symbol with lower index is defined as derivative with respect to the specified 
upper index quantity and vice versa. Representation (pSj) of the interaction term leads 
immediately to 

-iT'i,{6/5vfZ[v] + Tr^CZlv + C^] + Va^iv] = 0- (27) 

Finally introducing the generating functional W of connected Green's functions, defined 
through Z = exp(i>V), we arrive at 

i6/5r]rW[r]] = -Gfr], - G^^ Tr^C^^ exp(i5W^M) , (28) 

where 

SW"^ [r]] = W[r7 + C^] - yV[r]] . (29) 

This is our central result, on which approximation schemes are based. Equation (|28| ) 
determines directly the response to external fields. Further derivatives allow access to 
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linear response and correlation functions, and so on. A straightforward procedure to 
obtain systematic approximations is to truncate the expansion 

oo ^ 

5>V^W = ^ - Ci {^l^vr- {6l6r,r-W[r,]. (30) 

n=l 

The lowest order approximation, from now on referred to as mean-field approximation, 
keeps only the first term (n = 1) in ([30|). The quantities = {5/5riY'W and G""^ = 
— {5/5viY{5/5rD^yV, defined for r]^ = {Ex,0), are determined by ( pSf ) and its derivative: 

r = -GfVb - Gf[TT\^exp{iC<f>% (31) 
G-" = Gf + GS^[-iTr\,%^exp(iCV)]G'^^ (32) 

Straightforward calculation shows that the first equation is simply the saddle-point con- 
dition. In the second one the self-energy S^^, which equals the expression in squared 
brackets, is merely dependent on the mean value 0", and therefore neglects fluctuation ef- 
fects as well. We add here that the self-energy follows from Dyson equation T^f^G'^'^ = 5^^ 

through r,, = ro, - s,,. 

The second order approximation keeps the first two terms (n = 1 and n = 2) in (|30D. 
Unfortunately the analogous equations do not form a closed system. The third derivative 
of W couples to the self-energy. Omitting this, we obtain the so called self-consistent 
Hartree appoximation: 

r = -Gl\ - GS^[Tr^4^exp(iCf0^ - (i/2)Cf , (33) 

^ ^ab ^ G-[-iTr\,%^exp(iCV - (i/2)C^G%^)] (34) 

Again the expression in squared brackets in the second equation is identified as the 
self-energy. In contrast to the mean-field equations, the above approximation takes into 
account gaussian fluctuations in a self-consistent manner. 



5 Results and discussion 

We focus now on the case, where the CDW is moving uniformly with velocity w in a 
constant field E^, i.e. {{ip)) = vt. This translates into r]^ = {£^,0) and 0° = (0, ((v^))) 
for the representation chosen in the preceding section. Furthermore, G"'^, which solely 
depends on differences of time and space arguments, has the structure ([l6|) . 

As it turns out that mean-field equations give no corrections to free motion, except 
for correlations, we consider now the self-consistent Hartree approximation. Inspection 
of these equations shows immediately that the parameter A merely appears as a non- 
relevant short-time cut-off. Therefore we display equations (|33|j3^) in the classical limit 
A ^ oo: 

v = E^ + An [dr e'^^") sm{vT) G^(0, r). (35) 
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The self-energies are 

S^(r-, t) = 47r 6{r) [ V{t) G^(0, t) - J dr V (r) G^(0, r)] , (36) 

j:^{r,t) = -4m6{r)V{t), (37) 

where 

M(t) =i[G^(0,0) -G^(0,t)], (38) 

V{t) = e-^^W cos(t;t). (39) 

The quantity M(t) equals the correlation function (1/2) {{[6ip{0,t) — 5^9(0,0)]^)), where 
Sip = ip — {{(p)). We note that the selfenergies are only frequency dependent, so wave- 
vector integration can be executed in Fourier representations of and . After some 
algebra, the Fourier transform of (|36|) reads 



(40) 
where 

Z« = + _ (41) 

is specific to c? = 3. By convention, Im^ > 0, i.e. the cut is on the positive real axis. In 
addition we find 

M{t) = |^(l-e--)(-ir^)(87rImZ,^)-^ (42) 

recall that = F^ — S^. Thus the two self-energy components are interwoven in a 
non-linear fashion. 

From now on we neglect the temperature, i.e. we set F^ = 0. In doing so we do not 
miss any qualitative features, because long-range order is not destructed by thermal noise 
for d > 2. Numerical simplifications in this case stem from the fact that M{t) and V{t) 
are now periodic functions in time, without any superposed short time transients. This 
is reflected by the ansatz = —i{2n)'^^^VnS{uj — nv), where summation runs over 
all integers. Note that the n = contribution induces time persistent correlations. The 
numbers Vn = {v/n) f^^J^-dtV(t) e^""^^ are the Fourier coefficents of V{t) and determine 
self-consistently M(t) = J2nj^o^ni^ — e^™*)(4ImZ^^)^^ It turns out that quantities at 
frequencies nv, n E Z, are sufficient to form a closed set of equations. In this case Eq. (^) 
reads S^^ = (l/2i) ^^V^{Z^^^^^ - Z^J. Given v, the values of the K, ^nv or Z^^ 
respectively, are calculated by numerical iteration. Then the force- velocity dependence 
(|35D follows from 

v + {l/2)Y,Vn\^Z^, = E,. (43) 

n 

Due to the square-root in the expression (^T]), a small inertia /i is irrelevant at not 
too high velocities. That case is easily treated by neglecting the self-energies in a first 
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iteration of the self-consistent equations. For our purposes the hmit /i — is appropriate 
and used further on. 

The result is shown in Figs. |l] and where velocities v down to 10^^ are considered. 
The non-linear conductivity v/E^ develops an upward curvature below ~ 0.5 or 
V ~ 0.1, respectively. As shown by subsequent numerical results, this marks the crossover 
between the high- and low-velocity regime. A sharp drop in velocity appears at E^ 0.25 
in Fig. |. 

Further information is given by the equal-time phase correlation function 

2M(r) = (([(5v9(r,0)- (5(^(0, 0)]2)). (44) 

It is connected to the Green's function via M(r) = i[G^(0,0) - G^(r,0)]. The 
long-range behaviour, 2M(r) \r\/Li, is determined by the length Li = 2/Vo. Hence 
phase coherence is absent in the driven CDW as well. In the high-velocity regime, where 
Li ^ (8w)^/^, this spatial dependence is valid except for short distances |r| < 1. Figure 
^ shows Li as a function of velocity in the considered range, where Li is typically one 
order of magnitude larger than the bare Lq. It increases again at smaller velocities in a 
power law fashion, Li oc -y""'^^. 

The spatial dependence of phase correlations for small velocities is depicted in Fig. ^. 
Clearly three regimes as a function of distance appear: 

l^l for |r| < 1, 

2M(r) ^ <( 21n|r| for 1 < |r| < Li, (45) 
\r\/Li forLi<|r|. 

The prefactor two of the logarithm, here determined within an error smaller than one 
percent, is identical to that in the static CDW as calculated in pO|, |21| by the replica 
method. 

At high velocities the asymptotic behaviour of the normalized equal-time velocity 
correlation function, 

C(r) = ((<50(r, 0)50(0, 0)))/^;^ (46) 

is given by C{r) ~ L2 exp(— f) sin(f)/f, if expressed as a function of the rescaled distance 
r = \r\/L2. The new length L2 — {v/2)~^^'^, determines the decay of the oscillating ve- 
locity correlations. Thus, in contrast to the phase itself, velocity coherence is maintained 
in the sliding CDW. 

Inspection of C{r) in Fig. |^ shows that the situation is more complex in the low- 
velocity regime. Diverging relative on-site fluctuations with an approximate power law 
C(0) oc v"^'"^^ are found. In the following, we define a length L2 as in the high- velocity 
regime. On shorter distances, a velocity independent exponential attenuation is exhib- 
ited, which changes on larger distances into an exponential-like decay, superposed by 
oscillations. The long-range part of C{r) at different velocities can be projected on a 
single curve by rescaling distance through f = \r\/L2. Surprisingly, L2 exhibits a similar 
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power law as Li as a function of v, although the prefactor is somewhat smaller. A rather 
good choice for the overall functional form is 

C(,)^|^(0)-P(-|-l/2) forM<L„ ^^^^ 
I exp(— f) sin(r)/f for L2 < \r\. 

Our length L2 corresponds to ^ as introduced in Eq. (1.5) in From the static 
distortions of the phase (1.20) in p|, we read off the correspondence of our Li to 
where the exponent = is given in (1.21). This is in contradiction to our scenario, 
which would give 77^ = 1, and also to known results for the pinned CDW. 

The frequency spectrum of velocity correlation shows narrow-band-noise due to 
time persistent correlations: The zero wave-vector component of the Fourier trans- 
formed C(r,t) = {{6(p{r,t)6(f (0,0)))/ v"^, determines the noise spectrum S{uj) = 
J d^rdtC{r,t) exp{iujt). We obtain S{uj) = i{uj /vYGqJL^Gq^^. This expression con- 
tains only peaks at integer multiples of -u, i.e. S{uj) = (27r)^ S'„ ^(0; — nv), re- 
flecting the periodicity of the sliding state. Its first harmonic Si is oc v'"^ at high, and 
approximately cx: v~^-^''' at low velocities. 

Now we relate the quantities characterizing the velocity correlations to the velocity 
coherence length L2. Thus C(0) oc L2 at high, and oc L2 at low velocities. Correspond- 
ingly, 5*1 oc L2 at high, and oc L2 at low velocities, which is in agreement with the 
behaviour of dynamic correlations (1.18) in j^. 

In closing this section, we discuss the time dependence of force correlations. The 
interaction term (^) contains the fiuctuation component (f^ quadratically (for A 00) 



and the force-force factor cos{ip'^ —ip'^) in which the center-of-mass component y?^ appears. 
It follows that its mean value is proportional to the second derivative of the vertex 
functional with respect to (p^ or the self-energy component Sij respectively. Actually, 
V{t) = {{cos[ip{0,t) — ip{0,0)])) coincides with the definition (^). Its renormalization 
from the high-velocity limit V{t) ~ cos{vt) is shown in Fig. |^, where the peak at t = 
sharpens considerably at low velocities. Except for this singular behaviour, the primary 
dependence at low velocities appears to be a scaling of the global amplitude oc 
Whether the width of the peak vanishes on further lowering v has to remain open. 



6 Conclusions 

We have considered the sliding CDW in presence of quenched disorder within phase 
dynamics. Our starting point is a field theoretical technique, which allows disorder aver- 
aging directly in the generating functional, thereby introducing an effective interaction. 
To obtain explicit results, we have developed a self-consistent approximation scheme. 
In three space dimensions, the sliding regime at average velocity v is characterized by 
linearly growing phase fluctuations with spatial distance. The scale is set by the diverg- 
ing phase correlation length Li oc v^^'^ for high, and Li oc v~^-^'^ for low velocities. For 
low velocities a range with logarithmic increase with distance develops, as known from 
the pinned CDW. Nevertheless, velocity coherence is maintained, and the sliding CDW 
is accompanied by narrow-band-noise. Velocity correlations show oscillating decay with 
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distance under an exponential-like envelope. The corresponding length L2 is oc -u^^/^ 
for high, and oc v'^'^"^ for low velocities. The high-velocity behaviour of correlations is 
known from perturbation theory. For low velocities the exponents of Li and L2 coincide 
within few percent. In addition the above behaviour of the velocity correlation is only 
the asymptotics for long distances > L2, whereas for distances < L2 an exponential at- 
tenuation on a velocity independent characteristic length appears. The amplitude of the 
first harmonic of the narrow-band-noise normalized to the mean velocity is proportional 
to the coherence volume L\. 

Our approximation has identified a new regime at low velocities with different scaling 
properties, although it is possibly not applicable in the critical regime close to threshold. 
Direct comparison with experiment is difficult because of the possibility of phase slip and 
finite size effects in real materials, but numerical simulations appear to give a similar 
picture. Finally we represented a distinct approach to systematic approximations, and 
proved equivalence to our approximation scheme. 

Acknowledgements: We thank J. Miillers for several helpful discussions, and 
P. Schmitteckert for computational advice. 



12 



Appendix 

In this paper the pivotal element for generating systematic approximation is a functional 
differential equation for W. Evaluation proceeds in assuming that interaction induced 
fluctuations SW^ are restricted to a finite order Taylor expansion. First order expansion 
leads straightforwardly to a mean-field description. Whereas to second order an addi- 
tional approximation has to be executed, namely omitting third order derivatives of W, 
to close the equations. 

Recently Miillers und Schmid applied functional Legendre transform techniques 
to improve analysis of vortex dynamics in disordered type-II superconductors. Therein 
a Legendre transform is performed for both an external field t]^ and a quadratic source 
K^f^. Of course the essential approximation lies in the fact that the action functional 
is calculated to first order perturbation theory. In the following we show equivalence 
of the above to our self-consistent Hartree approximation. Starting from the generating 
functional 

W[v,K] = -iln^ DM exp(i5[^] + ir]a^^ + {i/2)^'^K ^,^') , (48) 

new variables 0" and G""^ are defined via 

{5/5riYW = 0^ {6/6Ky''W = (l/2)(iG"^ + 0"0^). (49) 

At the end of the calculation one sets K^f^ = to extract physical information. The 
desired action functional 

A[<P, G] = Wh K] - r/,r - {mrK,,<P' - {l/2)K,,G'^ (50) 

is the Legendre transform of (^8|) , where rj^ and K^^^ are to be eliminated using (49). 
The derivatives obey the inverse relations 

(5/50)a^ = -Va - K,,4>\ {S/SG),,A = -m)Kab- (51) 

In order to proceed the action = 5o[0] + is now composed of a free part 

Sq = (1/2)0°T°^0* and an interaction Si = — i Tr"^ exp(i^^(y9'*). It has been shown that 

A[<P, G] = S[<P] - In r°,G- - (i/2) (5/ - G'^\5/5<P\{5/5<P),S[<P] ) 

-iln y"DMexp((i/2)^«r,,/ + i52[0,¥.])|^^^, (52) 

where the functional integral is normalized to its free part, and the new interaction 
term iS2[0, v^] equals iSi[0 + ip\ minus its Taylor expansion up to second order in <^°. 
The index 2PI restricts a diagrammatic perturbation theory to two particle irreducible 
contributions. As usual, L^^G'"^ = 5^. Finally expanding the logarithm in (0) up to 
first order in ^2, the following result is found: 

^[0, G] = {l/2)rrl,<p' - (1/2)5," In r°,G- - (i/2)(C - r°,G"") 

-iTr^exp(iCV - mCG'^'Cb)- (53) 
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Performing the derivatives (jST]) of this quantity leads to the conditions ( ^3]) and (|3^) . 

To compare the two approaches, one may argue that both determine minimal self- 
consistent approximations including fluctuations. Whereas the above Legendre transform 
is truncated at first order in the interaction, our method restricts explicitly interaction 
induced fluctuations to gaussian ones. 
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Fig. |l]: Non-linear conductivity v/E^ vs. electric field in (a) self-consistent Hartree 
approximation, (b) Hartree approximation using the free Green's function Gq*, and (c) 
first order perturbation theory. 




Fig. ^ Average velocity v vs. electric field in self-consistent Hartree approximation. 
Note the logarithmic scale on the vertical axis. 
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Fig. ^: Correlation length Li vs. average velocity v in self-consistent Hartree approxima- 
tion. Asymptotic power laws are Li oc v"^-^"^ for low, and Li oc v^'^ for high velocities. 




Fig. ^: Static phase correlation 2M(r) = {{[5ip{r,Q) —5^9(0,0)]^)) vs. distance \r\ in 
self-consistent Hartree approximation. The velocities v are displayed in the inset. The 
middle part equals asymptotically 21n|r|. 
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Fig. 1^: Absolute value of the static velocity correlation C{r) = {{S(p{r,0)Sip{0,0))) /v"^ 
vs. distance |r| in self-consistent Hartree approximation. The velocities v are displayed 
in the inset. Sharp drops indicate a change of sign. 




Fig. ^ Absolute value of the renormalized force correlation V{t) vs. time over one period 
27r/v in self-consistent Hartree approximation. The high velocity limit is V{t) ~ cos{vt). 
The velocities v are displayed in the inset. Sharp drops indicate a change of sign. 
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